Chemical and Mechanical Signaling in Epithelial Spreading 



ON 

o 
o 

> 

o 

NO 

(N 



43 ' 

Oh! 
I 

o : 



CZ3 1 

o ■ 
>v 

43 : 

Oh. 



> 

oo 
o 
in 



On 
O 



X 



L. M. Pismen 1 and S. Y. Shvartsman 2 

1 Department of Chemical Engineering and Minerva Center for Nonlinear Physics of Complex Systems, 
Technion - Israel Institute of Technology, Haifa, Israel 
2 Department of Chemical Engineering and Lewis-Sigler Institute for Integrative Genomics, 
Princeton University, Princeton NJ, USA 

We propose a minimal mathematical model to explain long-range coordination of dynamics of multiple cells 
in epithelial spreading, which may be induced, under different conditions, by a chemical signal, or mechanical 
stress, or both. The model is based on chemo-mechanical interactions including a chemical effect of stress and 
concentration-dependent traction. The results, showing kinase concentration distribution and cell displacement, 
allow us to reproduce two activation waves on different time scales observed in the experiment, and distinguish 
between distinct dynamical patterns observed under conditions of injury or unconstraining. 



Regulated movement and deformation of epithelial layers, 
two-dimensional sheets of cells that tightly adhere to each 
other, is a key process in adult and developing tissues. Dur- 
ing embryogenesis, spreading and folding of epithelia plays a 
central role in gastrulation, an event that initiates the forma- 
tion of three-dimensional structures of tissues and organs. In 
adult tissues, epithelial migration is important for wound heal- 
ing. Mechanisms of collective cell migrations are a subject of 
intense research [1]. It is generally believed that epithelial 
spreading involves polarization of the cell cohort propagating 
from "leading" cells to the "followers". The process is medi- 
ated by the intracell actin machinery and proteins of adherens 
junctions, and regulated by complicated chemical pathways, 
which still remain largely unexplored. 

Early models of epidermal wound healing J2l, adopted also 
in a number of later publications, were based on reaction- 
diffusion equations describing cell motion and proliferation in 
response to a diffusive chemical signal secreted by the wound. 
This approach turned out to be inadequate for the descrip- 
tion of spreading epithelia, which has to include the effect 
of both chemical and mechanical factors on collective cell 
motion. Recent in vitro studies gave evidence of the impor- 
tance of chemo-mechanical interactions, whereby directional 
cell migration is maintained by a positive feedback loop be- 
tween cell deformation and intracellular signal transduction. 
Matsubayashi et al 13|] proposed that collective cell migration 
during wound healing is caused by an activation (phosphory- 
lation) wave of Mitogen Activated Protein Kinase (MAPK). 
The experiments of Nicolic et al J3l using canine kidney cells 
followed up on this by tracing both spreading dynamics and 
distribution of activated kinase levels under different bound- 
ary conditions. 

Wounding generates both free space for spreading cells, and 
an active species emitted as a result of injury that initiates 
the spreading process. Spreading also occurs, however, al- 
beit at a somewhat slower speed, in the absence of injury, be- 
ing initiated by a purely mechanical action following uncon- 
straining of the epithelial layer. In the case of spreading initi- 
ated by injury, two waves of MAPK activation were observed 
(Fig. [TJ). The first fast wave is excited by the signal emitted 
at the wounded edge; consequently, it rapidly decays, and is 
followed by a slow wave gradually spreading into the interior 



of the epithelial layer. In the absence of injury, only the slow 
wave is present; its form and spreading speed are similar to 
those in injured epithelium, indicating that the same mechan- 
ical signal is likely to be operating in both cases. The experi- 
ment confirmed as well the crucial role of MAPK, proved by 
suppression of spreading in the presence of a MAPK inhibitor 

m. 

Very recently, Trepat et al 15]] reported direct measurements 
of cell traction and stress in a spreading epithelial layer (in 
the absence of injury), which proved that traction is not re- 
stricted to a proximity of the leading edge, as has been com- 
monly believed, but remains substantial deep inside the layer. 
This implies an effective mechanism of propagation of active 
traction within the tissue. In this Letter, we propose that di- 
rectional collective cell migration is maintained by a positive 
feedback loop between cell deformation and intracellular sig- 
nal transduction. Thus, mechanical and/or chemical signals at 
the wound edge can coordinate shape changes and movements 
of multiple cells deep within the epithelium. 

Our aim is to find out and test a minimal phenomenological 
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FIG. 1 : Distance from the wound edge that a MAPK wave reaches at 
different time points after injury (scratching, solid line or membrane 
peel-off, dashed line) and unconstraining without injury (dotted line). 
In the latter case, only one wave is present. Insets: schematic repre- 
sentation of the MAPK waves reaching different distances at differ- 
ent times [4fl. 
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model compatible with available observations. While chemi- 
cal pathways and chemo-mechanical interactions are still un- 
known (and are certainly far more complex than the rough 
model to be described), there are few well-pronounced exper- 
imental features that need to be explained: (a) faster advance 
in the presence of injury, especially at the leading edge [4j ; (b) 
activation of MAPK in the absence of injury, and similarity of 
the emerging MAPK wave with the second wave in the pres- 
ence of injury [4]; (c) strong internal traction, large internal 
stress |Ht] and weak decay of the speed of advance with depth 
J3] in uninjured epithelia. 

Our main task is to discern a minimal set of chemo- 
mechanical interactions that would allow one to explain the 
above features. We will concentrate therefore on estimating 
forces and displacement velocities averaged in the direction 
parallel to the leading edge, and abstract at this stage from a 
strong random component of the cell motion (J, |6[] causing 
stress and velocity fluctuations, as well as transverse instabil- 
ities which naturally arise when positive feedback loops are 
operative flUl- 

Since the motion is strongly overdamped, the appropriate 
equation governing the displacement velocities v.- L should be 
based on Aristotelean rather than Newtonian dynamics, and 
has the form similar to the Darcy law: 



= Kfi 



fi = dj<7 



(1) 



Here fi is the force vector, which can be defined in a usual 
way through the stress tensor aif, K is the mobility coeffi- 
cient, or inverse friction factor, and dj = d/dxj is the spatial 
derivative; summation over repeated indices is presumed. Al- 
though the mobility should be generally anisotropic in a po- 
larizable medium [9], this is irrelevant in the current context 
of the study of averaged motion where the only non-vanishing 
components are those in the direction x normal to the edge of 
the tissue, / = f x and a — a xx . 

The forces acting in the epithelial tissue are at least of three 
kinds. First, is an external "wetting" force f^ w ' applied at the 
leading edge, responsible for a tendency of an unconstrained 
layer to spread to available space. The other passive force is 
due to the tissue elasticity. The simplest model is the Hooke's 
law /( e ) = E d x u where u is the displacement along the x 
axis and E is the elastic constant. Finally, there is an active 
traction force arising within the cells. 

MAPK activation, caused by either chemical or mechanical 
signaling, plays a crucial role in cell spreading, as indicated by 
suppression of spreading by MAPK inhibition [4]. Although 
the mechanism of MAPK action remains so far unclear, it is 
reasonable to assume that it enables active traction, so that the 
traction force can be assumed to be proportional to the level of 
MAPK activation m, = Fm. The experiment suggests 
that MAPK should be activated both chemically and mechan- 
ically, but mechanical activation is likely to be delayed, to ac- 
count for the observed retardation of the MAPK wave in the 
absence of injury (Fig. [TJ. We postulate therefore that MAPK 
is activated directly by the signal emitted by the injury at the 
leading edge, but its mechanical activation occurs through an 



intermediate chemical species, which is produced by straining 
stress. 

Since we strive to keep the model as simple as possible, we 
will assume linear kinetic relationships whenever it is pos- 
sible. One should, however, be cautious to avoid runaway 
in the model containing a positive feedback loop. Such a 
loop is created here by the stress arising due to the motion 
driven by MAPK-activated traction, which, in turn, activates 
more MAPK. The runaway instability is prevented by assum- 
ing the dependence of the intermediate species production rate 
on stress to be saturable. Thus, we assume the following ki- 
netic model governing the concentrations of the signal, s, and 
intermediate species, c, and the level of MAPK activation, to: 



Dd 2 x s - s, 
a<jH(a) 
1 + bo ~ 

s + c — km 



qc, 



(2) 
(3) 
(4) 



where H(a) is the Heaviside step function. The only diffusive 
species in this scheme is the signal, with the diffusivity D; the 
dot denotes the derivative with respect to the time variable, 
which is made dimensionless by using the inverse signal de- 
cay constant l/k s as the time scale; other decay constants are 
actually ratios, q — k.Jk s and k — k m /k s . The MAPK acti- 
vation constants in Eq. © are eliminated by adopting suitable 
relative concentration scales. 

Another problem to be addressed in model building is that 
of spreading of stress within the tissue. Recent experiments 
ifioll suggest that stress does not propagate inside cells or tis- 
sues instantaneously or with the speed of sound, but in a de- 
layed fashion, showing both appreciable delay and attenuation 
with distance, which may be compatible with the poroelastic 
model [11]. Lacking reliable data, we stopped at the elas- 
tic model (Q]) enhanced by advective propagation of traction 
against the direction of motion ||5|]. Accordingly, the contin- 
uum version of the mechanical model ([1]), collecting the elas- 
tic force f^ and traction f^ a \ reads 



"X( m — oid x m) + 7 d x u, 



(5) 



where \ — KF, 7 = KE, and a is the fraction of the cell 
traction transmitted backwards into the tissue, and the trac- 
tion is assumed to be negatively directed. The suitably scaled 
stress is computed by integrating the r.h.s. of Eq. @: 



<T = f 



(w) 



X 



mdx — am I — ^{dxU — 1). (6) 



The spatially discretized one-dimensional discrete version of 
Eqs. (0, © to be solved below is 

u„ = -Xn[m n - a(m n+ i - m„)] 

+ j{u n +l - 2u„ + Un-l), (7) 

<?n = f [w) + X m i ~ al7ln + 1 

- j(u n+ i -u n - 1). (8) 



where u n , n = 0, . . . , N is specified as the position of the 
front of nth cell, and the effective mobility Xn is allowed to 
be position-dependent. 

The first stage, present in a wounded tissue only, is chem- 
ical activation by a pulse emitted at the leading edge and dif- 
fusing into the tissue. The diffusion is very fast on the time 
scale of cell motion, reaching the depth of dozens of cells in 
a matter of minutes (see Fig. [TJ, and therefore can be treated 
without taking any reshaping of the tissue or its cellular struc- 
ture into consideration. In this approximation, the solution of 
Eq. (fJJ in the half-plane x > for a semi-infinite tissue with 
the instantaneous pulse at the leading edge s(x, 0) = 6(x) as 
the initial condition and the no-flux boundary condition at the 
leading edge x = is 



s(x, t) 



I 



exp 



x 2 \ 
1 ADt) ' 



(9) 



The solution is normalized by the total injected signal, so that 
J„ s(x, t) dx = e~*. The solution for the level of MAPK ac- 
tivation, neglecting at this stage the mechanical contribution, 
can be also expressed analytically: 



m(x, t) 



-kt 



2y/D(l-k) 




Xy/{l-k)/D 



, (10) 



where erf(...) is the error function. The analytical expressions 
(O, ( [Tol l are applicable in the Lagrangian frame with the co- 
ordinate x referring to their initial positions. Both the signal 
level and the level of MAPK activation pass through a maxi- 
mum with time at any location. Another possible assumption 
of a constant production rate above a certain threshold restricts 
MAPK to a finite band with a width first growing and then 
shrinking with time, and dynamics of its concentration within 
this band is qualitatively similar to the above. 

As follows from from experiment [4] and integrating 
Eq. (0 using m(x, t) given by Eq. (ITOb with the decay con- 
stants ratio k « 0.1 that fits the dynamics in Fig. \T\ the ad- 
vance at the purely chemical stage (ending with the decay of 
both the signal species and signal-activated MAPK) is limited 
to a few cell lengths. This feable advantage cannot account for 
the observed substantial difference between overall long-time 
cell drifts with and without injury, which is most pronounced 
close to the scratched edge. To emulate this difference, we had 
to modify the model by assuming that traction is enhanced 
in the regions reached by the signal, e.g. due to accumula- 
tion of an additional very slowly decaying promoter, so that 
X = X0O-+X1S) with 



- 1 „-x/Vd 



S(x,t) = / s(x,t)dt, lim S(x,t) = 




2VD 



(11) 



FIG. 2: The cell displacement Au n (above), extension On — Un 
u n -i — 1 (middle) and stress (below) in the absence of injury in 
50-cell chain. The distances are measured in the units of cell size 
The parameters: D = 400 A 2 , 7 = 5 A~\ k = 0.1, fci = 2, q 
0.05, x = 0-5 A a = 0.01, 6 = 1, a = 0.2, = 10" 3 . 
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FIG. 3: The map of the level of MAPK activation in spreading fol- 
lowing injury(xi = 5, other parameters as in Fig.[2] 
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FIG. 4: Coomparison of cell displacement at the end of computation 
runs with (upper curve) and without (lower curve) injury. 



Another modification of the basic model we found to be 
necessary is acceleration of the MAPK decay under the in- 
fluence of stress; otherwise excessive drift is observed in the 
back where large stress accumulates in our computations, 
as well as in the experiment [5]. Thus, we set in Eq. 
k = fco(l + k%a). Under this assumption, the model con- 
tains actually a feedforward loop connecting the stress and 
the level of MAPK activation: the former activates the latter 
through the intermediate species and suppresses it directly. 

We integrated Eq. (O together with Eq. ^ and the dis- 
cretized Eqs. ©-(HI for a chain of 50 cells with the contin- 
uation boundary conditions at the back end imitating a semi- 
infinite layer. The computations were carried out with the help 
of the Mathematica ODE solver. In the absence of injury, the 
solution of Eqs. (O-© with zero boundary conditions at the 



leading edge and zero initial conditions is trivial, and a weak 
"wetting" force f( w ' applied at the leading edge is necessary 
to initiate motion, which is consequently accelerated by the 
positive feedback loop. An example of computation starting 
from an unperturbed tissue is given in Fig. [2] The stress is ac- 
cumulated in the back, in accordance with the data of Trepat 
et al (5J. The extension is also maximal at the back, which 
may explain that it is there that most cell divisions occur Jg], 
even though divisions, a relatively rare event, are not incorpo- 
rated in this model. The advance approaches a quasistationary 
speed in the model semi-infinite tissue, though, of course, it 
would stop in a finite tissue with suppressed cell divisions. 

With an injury, the stress-dependent activation of MAPK 
follows the signal-induced stage. In accordance with the ex- 
periment (4J], the the level of MAPK activation, after reach- 
ing a maximum, decreases before going up again, as seen in 
Fig. [5] The difference between the drift of front and back 
cells is here, in accordance with the experiment, far more pro- 
nounced than in the absence of injury. The displacements 
at the end of computation runs with identical parameters are 
compared in Fig. [4] 

We anticipate that the above minimal model built under 
conditions of uncertainty will encourage further quantitative 
studies on chemo-mechanical interactions in spreading tis- 
sues, and will, in its turn, be adjusted and modified as more 
data become available. 
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